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Abstract 

A two-dimensional Pauli Hamiltonian describing the interaction of a 
neutral spin- 1/2 particle with a magnetic field having axial and second 
order symmetries, is considered. After separation of variables, the one- 
dimensional matrix Hamiltonian is analyzed from the point of view of 
supersymmetric quantum mechanics. Attention is paid to the discrete 
symmetries of the Hamiltonian and also to the Hamiltonian hierarchies 
generated by intertwining operators. The spectrum is studied by means 
of the associated matrix shape-invariance. The relation between the in- 
tertwining operators and the second order symmetries is established and 
the full set of ladder operators that complete the dynamical algebra is 
constructed. 
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1 Introduction 

In this work we will study a Pauli Hamiltonian describing the interaction of 
a neutral spin- 1/2 particle interacting with a magnetic field generated by an 
electric current-carrying straight wire. This system was introduced in Ref. pQ, 
where it was analyzed in the momentum space. Here, we will carry out a system- 
atic study in the configuration space based on the techniques of supersymmetric 
(SUSY) quantum mechanics [21 El ? sometimes referred to as the factorization 
method 4 , a technique that has been already used to study spin-1/2 Pauli 
equations (for different approaches, see for example 5 ). 

To motivate the specific form of the Pauli Hamiltonian for the present work, 
we will start with the formulation of its symmetry properties. We will see 
later that this additional information is also closely related to the factorization 
method. Thus, let us consider the Pauli Hamiltonian in the three-dimensional 
space 

P 2 

W 3 = ^ + /ia-B(x)+V(x) (!) 
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where B(x) is a magnetic field, p, is the magnetic moment of the particle, V(x) 
is a scalar potential, and a = (a\, a 2 , 03), <Jj being the Pauli matrices, that is, 



1 \ /0 1\ / -i \ ( 1 



a ° = 1 o 1 J ai = { 1 J ° 2 = ^ i J CT3 - v - 1 

' (2) 

We take the usual notation for the cartesian coordinates x = (x-y, x 2 , X3) and 
the momentum operators pk = —ihd/dxk = dxk- Now, we will impose some 
symmetries on this Hamiltonian in order to obtain the precise example we want 
to consider in this work. 



1.1 Rotational-translational symmetry along the a;3-axis 

We assume that the current-carrying wire is placed on the a^-axis, and we look 
for the systems Q allowing for first order symmetries of the form 

= p 3 + o(x) J 3 =X!p 2 - x 2 Pi + 6(x) = L 3 + 6(x) (3) 

where a(x) and 6(x) are Hermitian matrix-valued functions to be determined. 
By imposing that J3 and P3 are symmetries of H3, that is, [^3, H3] = [P3, H3] = 
0, we are led to the following explicit expressions (up to an equivalence): 

J 3 = -ihdg + I3u3 V 3 = -ihd X3 (4) 
B(x) = (f(p) cos[2/5(0 - 9 )}J(p) sin[2/3(0 - 9 )],g(p)) (5) 

y(x) = v{ P ) (6) 

where {p,9,X3) are cylindrical coordinates, the functions f(p),g(p),V(p) are 
arbitrary, and ft, 9o are free parameters. Therefore, we can decouple the problem 
as follows: a free motion in the a^-axis and another motion in the plane (ari, X2). 
In the sequel, we will restrict ourselves to the analysis of the system in the plane, 
with Hamiltonian 7Y 2 = Ti.3 — p\/{2m). 



1.2 Parabolic symmetry 

Next, we will also assume that there exists a second order symmetry of parabolic 
type, that is, with a leading second order term associated to parabolic coordi- 
nates [Hj, 

Si = L 3 pi + pi L 3 + A(x) • p + p ■ A(x) + N(x) (7) 

where A(x) and N(x) are matrix-valued vector and scalar Hermitian functions, 
respectively. If the magnetic term determined by f{p) is present, we arrive, up 
to an equivalence, to the same model already reported in £Q: 

B(x) = {x 2 /p\ -xi/p\ 0) F(x) = (8) 
J 3 = -ihdg + ^a 3 = L 3 + ^a 3 (9) 
Si=J 3Pl + Pl J 3 -px 2 <T-B(x). (10) 
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The commutator of J 3 and S\ gives another second order symmetry S 2 . 

S 2 =J 3 P2+P2J3+fJ-X l( T-B(x). (11) 

All these symmetries, together with the Hamiltonian H 2 , close the following 
quadratic algebra: 

[J 3 ,S 1 ]=iS 2 [J 3 ,S 2 ] = -iS 1 [S 1 ,S 2 ] = -4iH 2 J 3 . (12) 

As it was mentioned in , the symmetries S\ , £2 are similar to the components 
of the Laplace-Runge-Lenz vector for the Coulomb potential. 

Previous works have considered the quantum mechanical problem that we 
are dealing in this paper, mainly in the momentum space \7\, but also paying 
attention to some partial aspects in the configuration space ^8 . In the present 
work we want to address it from a self-contained, systematic and more complete 
point of view, including important properties not considered before. Basically, 
we will use the SUSY quantum mechanics approach with special emphasis on 
the shape-invariance of the model. 

This paper is organized as follows. In section 2 we will perform the sepa- 
ration of variables and we will obtain the discrete symmetries that are basic in 
the development of the following sections. Next, in section 3 we analyze the fac- 
torization and the shape-invariance properties of the radial equation. We study 
the ground states by means of this factorization in section 4, and the excited 
states in section 5. In section 6, we will study the relationship between the 
second-order symmetries and the intertwining operators entering the factoriza- 
tion. In section 7 we build the matrix ladder operators connecting eigenstates 
with different energies of the same Hamiltonian. Finally, in section 8 we con- 
clude this work with some remarks, stressing the most original results obtained 
in this paper. 



2 Separation of variables and discrete symme- 
tries 

Let us consider again the two-dimensional matrix Hamiltonian TL 2 obtained 
from ifT}. where we have excluded the free part along the X3-axis, and where 
the interaction is given by the magnetic field JSJl, thus, sharing the symmetries 
©~C3- I n order to simplify the expressions, we take 

h = 2m = 1 x^fixi y — fxx 2 r — fj, p (13) 
so that we will work with the following form of the Hamiltonian 

H=H 2 /fi 2 = -{dl + d 2 y ) + ^o l ~^a 2 (14) 
where the second order symmetries are 

j = 1,2. (15) 



3 



2.1 Separation of variables 

Since this Hamiltonian commutes with the operator J% given in @ , we can look 
for their common eigenfunctions 

The rotational symmetry can be used to separate variables in polar coordinates 
(r, 9) so that the Jz -eigenfunctions in (|TT)|l are given by 

9\(r,0) = Y x (6)F x (r) (17) 

with 



MO) = [ q e Wme ) Fx(r) = { j . (18) 

Now, taking into account the polar expression for the Laplacian A = + 
l/rd r + 1/r 2 dg, the eigenfunction equation for 7i takes the form 

- F' x '{r) - - F x (r) + ^^F x {r) - A <j 3 F x (r) - - a 2 F x (r) = EF x {r). (19) 

In order to eliminate the first order derivative term, we make the replacement 

F x (r) = r-V 2 $ A ( r ) (20) 
so that we finally have a one-dimensional matrix Schrodinger-like equation 

Hx$x(r) = (--^ + A " ~ 2 Af73 - - j $ A (r) =E$ x (r) (21) 
|_ ar z r z r j 

where we have introduced the notation Tlx for the radial part of the initial 
Hamiltonian TL. Notice that this equation is quite similar to the radial Schrodinger 
equation for a charged particle in a Coulomb potential 

H c t tt>(r) = -4>"{r) + i-Zi 0( r ) _ I ( r ) = E c 0(r) (22) 

where ^ is the orbital momentum. The discrete spectrum is given by the well 
known formula 

EC = - ^ + l n + iy B = ' 1 '- < 23 > 

2.2 Discrete symmetries 

In the following we will describe some discrete symmetries of the matrix equation 

GO- 
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2.2.1 Conjugation. 

Let us consider the antilinear operator 

C = <t 3 K (24) 

where JC is the complex conjugation operator, JC^(r) — $*(r). Then, it is 
immediate to check that this is a symmetry of TL\: 

H\C=CHx. (25) 

The eigenf unctions of C, up to a global phase factor, are of the form: 

(^) 



i 4>2 (r) 

where <fii{r) and (fo(r) are real functions. Hence, from now on, we will choose 
the eigenf unctions of the matrix equation Q21(l in the form 1)26(1. 

2.2.2 Reflection in A. 

It is very easy to see from l|21|) that 

V2Hx=H-xa 2 . (27) 
This means that the eigenfunction problem for H\ is equivalent to that of TL-\: 

H X <t> X = E<S> x W-A*-A = £*-A *-A = ^2*A- (28) 

Observe that these two discrete symmetries can be implemented also in the 
space of eigenfunctions of the two-dimensional equations fTH|l . If we want to 
compare these discrete symmetries with the ones of the Coulomb Hamiltonian 
HI of l]22p. we see that the conjugation is translated into the real character of 
the eigenfunctions, while the reflection property means that Tig — TL c _g + i- 



3 Factorizations and shape-invariance 

In this section we will investigate the factorization and the supersymmetrical 
properties of (|21|l . the radial part of Pauli equation. This one-dimensional 
2x2 matrix problem is interesting by itself, for an arbitrary value of A, but 
especially for the two-dimensional physical model (fH|l with A = 1/2 + to, to e Z. 
In particular, we will see that this matrix model obeys the shape-invariance 
properties. This is very interesting because, up to now, the shape-invariance 
was used as a very elegant algebraic approach in the solution of one and two- 
dimensional scalar spectral problems [21 El EH (HQ , while in this section it will be 
used to determine the spectrum of a one-dimensional matrix example. 

Following closely the well known factorization of the Coulomb Hamiltonian 
l|22|) . we propose here the following ansatz for this matrix case: 

f d A + Ba 3 \ ( d A + Bcr 3 \ 

Hx=[ — + -+Da 2 --r + -+Da 2 +7 (29) 

\ar r J \ dr r J 
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where A, B, D and 7 are real constants to be determined. Indeed, we find two 
different solutions having the following form 

Kx = L-L+ + 7A = L+_i^a-i + 7A-1 7A = 4(A ~\ /2)2 (30) 

where we have used the following notation for the operators 

Lt= d + (X + 1/2) - (1/2)0, _ (l/2) 7 d + Wx{r) (31) 



being W\(r) the matrix superpotentials. 

Of course, equations (|29[) - l|31[) are valid for any value of A ^ —1/2, and 
therefore we can build a hierarchy of Hamiltonians {TLx+m\, m € Z, as follows: 

T-i-X+rn = L x+m L X+m + 7A+m = L X+ m -l L X+rn-1 + lX+m-l- (32) 

These relations allow to define an algebra of operators of the Hamiltonian hi- 
erarchy The operators L x+m act as intertwining operators between two 
consecutive Hamiltonians TLx+ m and Ti.x+ m +i- 

L \ +m ^x+m = 'Hx+ m +iLX +m Hx+ m L x+m = L x+m Hx+ m +i- (33) 

These relations imply that the eigenfunctions of Ti.\ +rn can be obtained from 
those of Ti-x+m+i with the same eigenvalue by acting with L x+m , while the op- 
erator L~^ +m makes this connection in the opposite way. This correspondence 
between square-integrable eigenfunctions is one-to-one, except when L^ +m an- 
nihilates some eigenfunctions. As usual, relations l|32l) and mean that the 
Hamiltonians of the hierarchy {H\+ m } are shape-invariant: intertwined Hamil- 
tonians have the same shape but different values of the parameter. 

Now, let us see the effect of the discrete symmetries C and 02 on the inter- 
twining operators. It is easy to check that 

CLt +m =Lt +m C (34) 

and 

From the commutation (|34fl we see that the interwining operators keep the form 
of the eigenfunctions <& in (|26[l . Relation IMol) will have consequences for the 
ground states, as we will see in the next section. 

The hierarchies obtained from values A and A + mo, where mo £ Z, are the 
same, so we can get a family of different hierarchies characterised, for instance, 
by the values A £ [0,1/2]. Although the initial two-dimensional model l|14f> 
is related to the hierarchy with half-integer values of A, we will analyze the 
properties for any A. Thus, we can distinguish the following three classes of 
hierarchies. 
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• A = 0. We call it the 'integer hierarchy' {H m }, m e Z. It includes the 
Hamiltonian H.q and would correspond to two-valued eigenfunctions in 
the context of the two-dimensional problem 114(1. The reflection in A is 
implemented in the hierarchy, a scheme of which can be seen in Figure ^ 

• A = 1/2. We refer to this case as the 'physical hierarchy' {H m +i/2}i m G 
Z, because it is associated to the physical spin-1/2 systems of Eq. (I14|) . 
Indeed, it is the only case considered in all the previous references [UI7IIH]; 
and for it the operators L_ t , 2 are not well defined, as can be seen from 

Therefore, there are no first-order intertwining operators connect- 
ing the Hamiltonians Ti-i/2 and 7i_!/ 2 - Remark that the hierarchy also 
includes the reflection in A. 

• 0<A<l/2. We will refer to this as the 'general hierarchy' {TL\+m}- It 
does not implement the reflection in A; in fact this reflection gives rise to 
the hierarchies with values —1/2 < A < 0. The associated two-dimensional 
eigenfunctions are multiple-valued. 



4 The ground states of the hierarchies 

It is known [21 that the ground states of SUSY-hierarchy Hamiltonians are the 
main elements when using the shape-invariance approach in the construction of 
the whole spectra and the eigenfunctions for one-dimensional scalar problems. 
In this section, the explicit expressions for the ground states of the matrix 
SUSY-hierachies described in the previous section will be found and analyzed. 

Let us start with the ground-state wave-functions $\ +m annihilated by the 
general intertwining operator L~^, m : 

Lt +m K +m = 0. (36) 

Then, according to lj3H)l and l|3*2*j l. this will be an eigenfunction of H\+ m with 
energy E° x+m = 7 A+m : 

T-L\+ni®°\ +m = £°+ TO $ A+m E \+ m = - ^x + m+ 1/2) 2 ' ^ 37 ^ 
Let us use the notation 1(26(1 for the ground state 



"o A+ ", (38) 



and let us make the substitution 

^« = ^W J ' = 1 ' 2 * = 2(A + m + l/2) (39) 

in equation l(36() . Then, we arrive to a modified Bcssel equation for the compo- 
nent fi{z) 

«V'(*)+M(*)-(* 2 + l)Pi(*) = (40) 
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with general solution 

i Pl {z)=c x h{z) + c 2 K l {z) (41) 

where Ii(z), K\{z) are the two linearly independent modified Bessel functions. 
After a simple calculation, we arrive to the general form of the ground state 
(|36|) for any integer m: 



with 



= c i K A+m(r) + c 2 I A+W (r) (42) 



KxWr) = '1»-tW (43) 



z/ ^ I 2(A+m+l/2) 



u +m (r)=r^+™[ 7l V^+V2)/ , (44) 



^ 2 (A+m+l/2) / 

Remark that as the asymptotic behaviour of the modified Bessel functions 
/o( r ) an d Ii( r ) m the limit r — > +oo is exponentially increasing, they can not 
lead to square- integrable functions. On the other side, in the same limit K$(r) 
and K\{r) decrease exponentially, while near the origin K$(r) m log(r) and 
K\(r) sa 1/r. Thus, the physical solution will be ^ +m (r) oc K A+m (r), which is 
a well behaved ground state of Tt\+ m with A + rn > (the ground state of Ho 
will be considered later). 

Let us consider now the ground states $° +m of H\+ m with A + m < 
0. They are annihilated by the operator L x+m _ 1 in the factorization (|30|l : 
^A+m-i^A+m = 0- Then, according to (|3*U)l and lO, this will be an eigenfunc- 
tion with energy E° x+m = E°_ x _ m = j-x- m : 

Wx+m^A+m = ^A+m^A+m ^A+m = ~~ _|_ m _ 1/2) 2 ' 

After a computation similar to the one carried out in the previous case, we get 
the following general solutions: 



*\+m(r) = CiT 




(46) 



As in the previous case, a study of the asymptotic behaviour shows that the 
square-integrable and finite ground states are obtained by taking C2 = in 
Ij46(l . whenever A + m < 0. The results for this case could be found directly 



8 



from the previous one with the help of the relation (|35fl between both types of 
intertwining operators through 02 • 

In summary, we have one reasonable ground state for each Hamiltonian 
7~L\+m'- if A + to > it is of the form &\_L_ m {r), while if A + m < it will 
take the form .^(r). There is only one Hamiltonian, Hq, having a doubly 
degenerated square-integrable solution with energy E® = — 1. But this case is 
rather special and we will comment on its lack of physical meaning later on. 

To end this section, let us discuss now the relationship between these ground- 
state solutions and the superpotential. First of all, let us recall that 

L+ +m K A+m (r) = L+ +m l x+m (r) = 0. (47) 

Let us introduce now the 2x2 matrix M\ +m (r), whose columns are K\ +m (r) 
and I> +m (r), in this order. Then, equations Q47[) can be expressed in a more 
compact form: 

L+ +m M x+m (r) = -M' x+m (r) + W x+m (r) M A+m (r) - 0. (48) 

Thus, after substituting 142(1. we get again the expression of for the matrix 
superpotential in terms of the solution matrix: 

W x+m (r) = M' x+m (r) M- X l m (r) = ( X + m +^-^ _ (4Q) 

Therefore, we have shown that this class of matrix solutions MA+mM gives rise 
in a nontrivial way to a Hermitian matrix superpotential. It is important to 
stress that there are not many explicit examples satisfying this condition 12 . 



5 Excited states 

In this section we will study the excited states of the system J2H- There are 
some general expressions valid for any hierarchy. In order to find the n-excited 
state ^(r) of a particular Hamiltonian TC\, A > 0, we can start with the ground 
state $° +n (r), n G Z + , of the Hamiltonian Tt\+ n in the same hierarchy. Then 

n(r) = L-L- +1 ---L- +n _ 1 ^ x+n (r) (50) 

with energy 

E " = ~4(A + n+ 1/2)2" ( 51 ) 
Two natural questions arise now: 

a) Are the excited states obtained in this way always square-integrable? 

b) Are these excited states the only bounded physical states for each Hamil- 
tonian of the hierarchy? 
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Here, in the context of the factorization method, we can answer to the first 
question: under the assumption A > 0, all the cigenfunctions build in the form 
(|5U|I are finite, square-integrable and vanishing at the endpoints (0 and +00), 
so they indeed describe excited bound states. The details are given in the 
Appendix. With respect to the second question, as in the two-dimensional 
scalar models ^U], it should be studied from the point of view of a general 
"oscillation theorem" for matrix Sturm-Liouville operators to guarantee that 
no additional excited states exists besides those obtained after applying the 
shape-invariance procedure (for a discussion on two-dimensional scalar models 
see JD]). Meanwhile, one can notice that each component of the constructed 
n-th excited states i|5U)l has exactly (n — 1) nodes. 

If A < 0, the excited states with similar properties are obtained from the 
other class of ground states: 

n(r) = Lti^-a ' ' ' L t-n *A-nW (52) 

with energy 

= - 4(A-n-l/2)* - (53) 
Now, we will consider some specific features of each hierarchy. 

5.1 'Integer hierarchy' (A = 0): H m , m e Z. 

A representation of this hierarchy can be seen in Figure It includes the 
Hamiltonian TIq which has the peculiarity that both procedures lead to two 
square-integrable excited states for each energy 



Therefore, in principle, each eigenvalue level is doubly degenerated. However, 
an important difference is that in each eigenfunction one of the components 
vanish at the origin, while the other one does not. This feature leads to some 
problems about the physical interpretation, that we will discuss now. 



The matrix Hamiltonian TIq in 121(1 can be diagonalized by means of the 
unitary matrix U — r^g^co + 

*o* = {-£ - t} * ~ = {-£ - 7} * d (55) 

where H d = UH Q W and $ d = Thus, we can write 

^° + ) (56) 
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nj n_j n_j n_j n_j n_j n_j 

rt_3 rt_2 rt—i /to rti rt-2 riz 



Figure 1: Energy level diagram for the 'integer hierarchy', including the infor- 
mation about the ground states and the operators L^. 

,2 -. 

where h± = —jjpz ± - are scalar Hamiltonians. In particular, the lowest energy 
(E = —1) normalizable solutions of the Schrodinger equation with Hamiltonian 
are as follows 



One can notice that these solutions, as well as the higher energy solutions, do 
not vanish at the origin, and therefore correspond to divergent mean values of 
the potential energy terms. A plot of some of these functions can be seen in 
Figure |21 By this reason, solutions |JSSJ, though being normalizable, do not 
belong to the physical sector of the model. 

Nevertheless, these non-physical solutions are quite useful from the point 
of view of SUSY quantum mechanics. Indeed, they are used to generate the 
intertwining operators to get a hierarchy TL m of shape-invariant non-diagonal 
Hamiltonians starting from TLq, which is diagonal. Besides, any TL m , ra^O has 
a sensible physical discrete spectrum with eigenfunctions vanishing at the origin. 
By extending the label m to a real number we get all the hierarchies TL\+ m - 

5.2 'Physical hierarchy' (A = 1/2): H m+1/2 , m G Z. 

The general scheme <|5U|) ~(I53 [1 is still valid here. But since in this hierarchy the 
operators L^,^ are not defined, in order to connect the excited states of Tii/2 
and H-i/2 we need a zero-order intertwining operator given by the reflection 




(57) 





(58) 



(23) 



$" /2 OC 0- 2 $ 



-1/2- 



(59) 
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Figure 2: Plot of the two components of the eigenfunctions (the first component 
in solid line, the second one in dashed line) for the ground and the first excited 
states of the Hamiltonians Hq (left) and H1/2 (right). 



It is interesting to remark that it is possible to include in a very natural way 
the operator 172 inside the set of intertwining operators L^ +m if we 'normalize' 
them by a factor: 

iA +m -(A + m + l/2)i± +m . (60) 

Then, the set {L^ +m }, m S Z, will act also as intertwining operators of the 
hierarchy 7i\ +m as in . but they are always well defined, and in particular 
i^ X j 2 = — <J 2) a s it should be. 

However, the expression of the Hamiltonians in terms of the normalized oper- 
ators are changed by these factors and it gives rise to the following 'commutation 
rules': 

Li +m -iLx +m -i - L- x+m L+ +m = -2(A + rn)H x+m . (61) 



5.3 'General hierarchy' H\ +m , m £ Z. 

If < A < 1/2 we will have two kinds of Hamiltonians in the hierarchy. The 
right-hand Hamiltonians 7i\+ m , m = 0,1,..., and the left-hand ones H\- m , 
m £ N. We have two special features: (a) The excited states of the right (left) 
Hamiltonians are only obtained from the right (left) ground states. However 
i^Li do not connect these two kinds of Hamiltonians, neither the reflection 
matrix 02 can be used to connect the two Hamiltonian sectors, (b) The spectrum 
of these two types of Hamiltonians are different. In this case the operators 
Lf +m , acting on physical states, do not always generate again physical states. 
A schematic diagram of this case can be seen in Figure 01 



6 Factorization operators and second order sym- 
metries 

As we have seen above, the action of the operators L^ +m on eigenfunctions of 
the Hamiltonian hierarchy change the label A + m in one unit, while keeping 
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A-4 rl\-3 rtx-2 rt\ ttA + l «A + 2 rt A + 3 



Figure 3: Energy level diagram for the 'general hierarchy', including the infor- 
mation about the ground states and the operators L x+m . 

the energy. We can write this property as follows 

Li +m ^\+ m {r) cc^x+m+i{r) i Hra $A+m+i(r) oc $> +m (r). (62) 

If we recall that the label A + m is for the J3 -eigenvalue, we can say that the 
operators L^ +m act essentially as lowering or raising operators for J%. However, 
when they act on the eigenstates of the Hamiltonians Ti.\ they preserve the 
energy eigenvalue. 

Now, if we go back to the second order symmetries 71,72 of the initial 
Hamiltonian Tt, and we introduce the operators T ± = ±i(71 =p iT<x)J1, they 
satisfy the commutation rules 

\Ji,T*\ = ±T ± [W.T*] = [T+,T-\ = -2HJ 3 . (63) 

This means that acting on the common eigenfunctions of Tt and J3 realises 
the same role played by L^ +m : change the J3 -eigenvalue in one unit while 
they leave that of Ti unaltered. Therefore, there must be a close relationship 
between both types of operators, as it was the case for the Coulomb problem 
|13| . Indeed, if we write the eigenfunctions ^\(r,8) given in 11711 H8|) in terms 
of $xW introduced in lj2U|) . 

9x{r,0) =r- 1 ' 2 Y x {9) $ A (r) (64) 

and we use the first commutation relation in Eq. I|63|l . we have 

T ± y x {r,6)<xy x ±i{r,0). (65) 

In order to prove this relationship explicitly, let us use first <|10f) - (|ll[l an d (|13|l - 
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I|15|) to express in polar coordinates: 



r~ = 



-ie 1 
2 



2 -3 9 + r ^ + -(a 3 +<7o) 



-2 —de + d r d e + -(a 3 + a ) - 



Then, if we take into account that 



3 r * A (r, 0) = r- 1 ' 2 Y x {6) + d r ^j * A (r) 



-e~ 
e l9 



tf A (r,0) = (-z)r- 1 / 2 F A (0)<7 2 $ A (r) 



(66) 

(67) 
(68) 
(69) 



and we insert these identities in l|66|) , we get the induced action of T ± on $ A (r) : 



T±$ A (r) = (A + ~ 



d (A ±1/2)- (1/2) a 3 1/2 
T dr r A±1/2 CT2 



$aW. (70) 



If we compare the above expression with the definition l|31|l of L^, we see that 
indeed, the action of the second order symmetries 7~ ± induced in the space of the 
eigenfunctions $>\+ m (r), essentially coincide with the action of the intertwining 
operators £ A+m , as we expected from the initial intuitive arguments. In fact, 

we get the normalized operators defined in l|60|) . and the commutation l|63|l 
of T induces in the space of functions $ A + m (r) the relation displayed in Q61JI 
between the operators £ A+m - 



7 Ladder operators 

In this section we will investigate the construction of ladder operators that will 
allow to link the excited states of a given Hamiltonian TL\. 

Let us consider the eigenvalue equation l|21|) for a given Hamiltonian 7i\ 
corresponding to the energy E™ and eigenfunction $ A (r) given in (|5(J|> — (|51|) . 
After multiplying by r 2 , we reorder the terms in the following way: 



dr 2 



rer 2 



4(A 



1/2) 



A<7 3 



*X(r) = -A 2 ^(r). 



(71) 

In a similar way to what we did in l|29|> - (|3U[l . we can factorize the operators 
1Z\ „ in two ways 



K\,n - Qt, n Q\,n + - <9 A ,„+1<2 A) „+1 + ^X,n+1 



(72) 
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where 



Qtn = 
Q\,n = P A+n 



dr 



^3 >,1 0~2 

2(A + n) 1 J 2 2(A + n+l/2) 



r ^ + 2(AT^- (A + n) -2 



<7 2 



2(A + n + l/2) 



" A '" = -4(AW- (A + n)2 + I 



and the operators T>\ 



are dilation operators defined by 
A + n + 1/2 / E\ 



A + n-l/2 



£1' 



(74) 
(75) 

(76) 



The operators Q^n+i ac * as intertwining of the diferential operators lZ\^ n and 

7^A,n+l : 

^A.n+l = ^A,n+1 Q\, n+ i Q \, n +l ^A,n+1 = ^A,n Q^n+f ( 77 ) 

As a consequence, the set of operators {Q\ n , n € Z + } will act as lowering and 
raising operators for the Hamiltonian H\, 

^A,n+l • V A ^ V A ^A.n+l ■ V A 



(78) 



From l|72[l we can compute the normalization of the eigenfunctions obtained in 
this way: 

||^ +1 || 2 = Ko-a, A , n+1 )||^|| 2 . (79) 

In particular for n — 0, we have Qq = 0. The operators Q^ n constitute a 
non-trivial generalization of those well known for the Coulomb problem . 

Finally, we can join the two families of 'ladder operators' found up to now, 
{^A+mi ano - {Qn }j to build a 'dynamical algebra' inside the hierarchy {H.\+ m }. 
However, only for the case A = 1/2 all these operators will connect exclusively 
physical states of the discrete spectrum. 



8 Remarks and conclusions 

In this paper we have studied a two-dimensional Pauli Hamiltonian, which has 
two independent integrals of motion, from the point of view of supersymmetric 
quantum mechanics. We have examined a series of properties that were not 
fully explored up to now: 

(i) We have considered some useful discrete symmetries since the very begin- 
ning. 

(ii) We also included in our study a family of Hamiltonian hierarchies labelled 
by the parameter AG [0, 1/2]. For all these hierarchies the spectral prob- 
lem was solved by means of the matrix shape-invariant method. 
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(iii) We have shown the relation between these intertwining operators and the 
second order symmetries. 

(iv) We also computed the ladder operators suitable to this matrix problem, 
which link the excited states of the same Hamiltonian. 

In summary, we have shown that the methods employed along of this paper 
(SUSY quantum mechanics) constitute a very useful tool for the investigation 
of Pauli matrix Hamiltonians, and deserve to be exploited for more general 
situations. 



Appendix 

In this appendix we will show that the excited physical states are well 

behaved near the origin and at infinity. We start with the expression (|50|l 
for the excited states, taking into account that the physical ground states are 
*A+n( r ) 00 K A+nM, where K\ +n (r) is given in |j33J|. For the sake of simplicity, 
we will use in the sequel the variable a — A + n, and therefore l|5UI) takes the 
form 

K-n(r) = L-_ n L-_ n+l ■ ■ ■ L-_ 2 L-_, S°(r) (80) 
where, according to l|43|) . the ground state is given by 

Using the new variable z, the operators adopt the form 



6 2a + 1 



d 2b+l-a 3 2a +1 

I (79 

dz 2z 2b +1 



(82) 



It is interesting to know the approximate form that have these operators near 
the origin (r, z ss 0) and at infinity (r, z ~ oo). Indeed, we have the following: 

it wo, ^-"^^(^ Ja) m 

if r.^oo, ^-^-^(^ If)- «») 

A.l Behaviour near the origin 

In this approximation, we have that 

®2-n~N-_ n N-_ n+1 ---N-_ 2 N-_ 1 $° a n = l,2,... (85) 

and using 

K (z) » log(z) K^z) » 1/z (86) 

K' (z) = -K^z) zK' 1 (z)=-zK (z)-K 1 (z) (87) 
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we can prove that 



2a- 1 \ f z a K x {z) 
2a + 1 J I iz a K (z) 



AT-_ 1 <«(2a+l) a 

^V.-2^a-i*a~^o+i; ^ Q [2a+l][2a-l] ) { iz^Koiz) ) 



Therefore 

rfo+i/2) 



^ a + ij I o r(a+3/2) I iz a ~ n+1 K (z) j 

\ r(a+3/2— "1 / \ ' 



(88) 



or 

r(A+ra+l/2) 



* S « 2"(2A+2n+l)^ r(A+n+3/2) • (89) 

\ u r(A+3/2) / V my/ 

From here, taking into account (|86[1 . it is easy to conclude that for any A > 
the states <£>™ are well behaved near the origin. 

A. 2 Behaviour at the infinity 

For big values of the variables z or r, we have that 

K-n « G"_ n G-_ n+1 • • • G-^G-^ $° (90) 
and taking into account that in this region we can approximate 

K u {z) « e" z (91) 

after a simple calculation we get 

_ (-<x 2 ) n r(q-n + l/2)r(2a+l) 

a ~™~ (2a + 1)™ r(o + l/2)r(2a-n + l) a 

or 

„ ^ {-<?2) n r(A + l/2)r(2A + 2n+l) 

A ~ (2A + 2n+l)« r(A + n + l/2)r(2A + n + l) A+n ' 1 ' 

Using (|91|l . we conclude that for any A > the states <fr™ are well behaved at 
infinity. 
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